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ABSTRACT 



Context. The presence of celestial companions means that any planet may be subject to three kinds of harmonic mechanical forcing: 
tides, precession/nutation, and libration. These forcings can generate flows in internal fluid layers, such as fluid cores and subsurface 
oceans, whose dynamics then significantly differ from solid body rotation. In particular, tides in non-synchronized bodies and libration 
in synchronized ones are known to be capable of exciting the so-called elliptical instability, i.e. a generic instability corresponding to 
the destabilization of two-dimensional flows with elliptical streamlines, leading to three-dimensional turbulence. 
Aims. We aim here at confirming the relevance of such an elliptical instability in terrestrial bodies by determining its growth rate, as 
well as its consequences on energy dissipation, on magnetic field induction, and on heat flux fluctuations on planetaiy scales. 
Methods. Previous studies and theoretical results for the elliptical instability are re-evaluated and extended to cope with an astro- 
physical context. In particular, generic analytical expressions of the elliptical instability growth rate are obtained using a local WKB 
approach, simultaneously considering for the first time (i) a local temperature gradient due to an imposed temperature contrast across 
the considered layer or to the presence of a volumic heat source and (ii) an imposed magnetic field along the rotation axis, coming 
from an external source. 

Results. The theoretical results are applied to the telluric planets and moons of the solar system as well as to three Super-Earths: 55 
CnC e, CoRoT-7b, and GJ 1214b. For the tide-driven elliptical instability in non- synchronized bodies, only the Early Earth core is 
shown to be clearly unstable. For the libration-driven elliptical instability in synchronized bodies, the core of lo is shown to be stable, 
contraiy to previously thoughts, whereas Europa, 55 CnC e, CoRoT-7b and GJ 1214b cores can be unstable. The subsurface ocean of 
Europa is slightly unstable. However, these present states do not preclude more unstable situations in the past. 

Key words. Hydrodynamics - Instabilities - Planets and satellites: interiors - Planets and satellites: dynamical evolution and stability 



1. Introduction 

The flows in fluid layers of planets and moons are of major inter- 
est because they imply first order consequences for their internal 
dynamics and orbital evolutions. Indeed, internal flows create 
torques on solid layers and induce energy dissipation. Moreover, 
internal flows are directly responsible for the generation of mag- 
netic fields, either by induction of an existing background mag- 
netic field or by excitation of a self-sustained dynamo. Finally, 
planetary heat fluxes are also directly linked to flows in fluid lay- 
ers, which can act as thermal blankets for stably stratified con- 
figurations, or as efficient heat flux conveyers in the case of con- 
vective flows. 

Planetary fluid layers are subject to body rotation, which 
implies that iner tial waves can propagate through them (e.g. 
iGreenspanI Il968h . Usually damped by viscosity, these waves 
can, however, be excited by longitudinal libration, preces- 
sion, and tides, which are harmonic mechanical forcings of 
azimuthal periodicity m - 0, 1, and 2, respectively. The 
fluid response to such forcings in ellipsoid s is a long-standing 
issue: see e.g. for longitudinal libration Aldridge & Toomre 
(|1969|) . [Noir et al . (2009), Calkins e t al. (20 1 0),"S aui-et et al.. 
J20ld) . IChan et al . (2011a) and Zh ang et al.l (1201 ll) . for lati- 

recession, IPoincare 



tudina l librati on, than et al ., (,201 Ih). for p rece 
(ll910l) . lBusse (1968). Cebr on et all fcoiOM . lKida & Nakazawa 
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("201^; "Kidal (wTlS and 'Zhan g et alj (l20Toh and for tides , 
Ogilvie & Lin (2004), Ogilvie & Lid (l2007h. iTilgnd (l2007h . 
Rie utord & Valdettard (120101) and lMorize et al.l (12010 ') In these 
studies, it has been shown that the dynamics of a fluid layer is 
completely modified when the forcing resonates with an iner- 
tial wave. In addition to these direct forcings, inertial waves can 
also form triadic resonances, leading to parametiic inertial in- 
stabilities. For insta nce, the s o-called s hear instabilitv can be ex- 
cited by precession (lKerswellLl993b;.Lorenzani & TiIgneiJ200lL 

l2003h . and the elliptical instability can be excited by tides in 

J ir — — 1 1 II — "f^ 

non- synchronized bodies (Malkus 1989; Rieutord 2000) and by 

librations in synchronized ones (Kerswell & Malkus 1998). 

The elliptical instability is a generic instability that affects 
any rotating fluid whose streamlines are elliptically deformed 
(see the review by Kerswell 2002). A fully three-dimensional 
turbulent flow is excited in the bulk as soon as (i) the ratio be- 
tween the ellipticity of the streamlines /3 and the square root of 
the Ekman number E (which represents the ratio between the 
viscous over the Coiiolis forces) is more than a critical value 
on the order of one and (ii) as soon as a difference in angular 
velocity exists between the mean rotation rate of the fluid and 
the elliptical distortion. In a planetary context, the ellipticity of 
streamlines is related to the gravitational deformation of all lay- 
ers of the considered body, coming from the static and periodic 
terms of the tidal potential, as well as from a potential frozen 
bulge. The differential rotation between the fluid and the ellipti- 
cal distortion can be oscillatory when caused by libration in syn- 
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Fig. 1. Sketch of the problem studied in this work. With the mean 
rotation rate of the fluid Q, we define the dimensionless orbital 
rotation rate 7(f) and the dimensionless spin rotation rate 1 + W(,). 
The phase lag between the tide and the gravitational potential of 
the host body is not relevant to our purpose. 



Table 1. List of the different astrophysical configurations that 
could lead to an elliptical instability (E.I) in a planetary fluid 
layer (liquid core, subsurface ocean) of a non synchronized or a 
synchronized celestial body. 



State 


Origin 

of Afi 


Origin 
of/3 


AQ 


E.I 


Non-sync. 


spin rotation 


D.L" 


^^spin ~ ^^orh 


TDEI 


Non-sync. 


impact 


S.B." 


spin-up process 


TDEI 


Sync. 


forced O.L* 


D.L 


2eQ.cos(2nt/T„ri,) 


LDEI 


Sync. 


forced P.L* 


S.B. 


eQ.co&{lTrtlT„h) 


LDEI 


Sync. 


free P.L 


S.B. 


eQ.cos{(i)f„et) 


LDEI 


Sync. 


any P.L 


S.B., D.T 


zonal flow^ 


TDEI 



yS is the ellipticity of the boundaries distortion, Afi is the differential 

rotation rate between the fluid and the elliptical deformation, il is the 
mean spin rate of the planet, e is the physical libration amplitude, 
Ti,rt and e are the orbital period and eccentricity, and ujfree is the free 
libration frequency. 

" D.T. stands for dynamic tides and S.B stands for static bulge. 
*■ O.L. and P.L stand respectively for optical and physical libration. 
^ Case equivalent to a non-synchronized case (cf. section l2.2.2t . 



chronized systems, or stationary in non-synchronized ones. The 
elliptical instability is then respectively refered to as libration- 
driven elliptical instability (LDEI) and tide-diiven elliptical in- 



stability (TDEI). TDEI and LDEI have al ready been suggested 

:al.iri997h 



as taking plac e respectively on Earth (e. g. lAldridge et ; 
and on lo (e.g. lKerswell & MalkuslI99a) . However, these previ- 
ous works do not consider some pl anetary particularities , so they 
need to be revisited. For instance, lAldridge et al.l (Il997b did not 
take the orbital rate of the Moon into account or the magnetic 
field of the Earth, thus neglecting the effect s of tide rotation and 
Joule dissipation on the growth of TDEI. iKerswell & Malkiia 
( Il998l) implicitly assumed that the tidal response of lo is com- 
pletely fluid, neglecting the rigidity of its mantle and overesti- 
mating the amplitude of librations and tidal deformations. Our 
purpose here is to extend previous results of the literature on 
TDEI and LDEI and to determine general formulas for quantify- 
ing the presence of the elliptical instability in terrestrial bodies, 
taking the relevant complexities present in natural systems into 
account. 

This paper is organized as follows. Section|2]presents the dif- 
ferent celestial forcings that could excite an elliptical instability, 
first focusing on tides in non-synchronized systems, and then on 
forced and free libration in synchronized ones. In section |3] we 
introduce our physical model and develop a local WKB analysis 
in all configurations, including the effects of viscosity, as well 
as the effects of an imposed magnetic field and a local temper- 
ature gradient. These theoretical results are used in section|4]to 
investigate the possible presence of any of TDEI and LDEI in 
telluric planets and moons of the solar system, as well as in two 
Super-Earths of extrasolar systems. The possible consequences 
of those instabilities are finally considered. 



2. From celestial mechanics to the excitation of an 
elliptical instability 

Figure [T] presents a sketch of the problem considered in this 
work. We consider a telluric body of rotation rate O.,^,,,, orbiting 
around an attractor (orbit in black dashed line) at the orbital rate 



Q.orb- This body has a radius R, a mass M, and a fluid layer in its 
interior between an external radius R2 and an internal radius R\, 
typically a liquid outer core. We suppose that this internal fluid 
layer is enclosed between an external elliptically deformed solid 
layer and a possible inner core, such as the outer liquid core of 
the Earth or the subsurface ocean of Europa. The elliptical de- 
formation can admit different origins. First, in the presence of an 
orbiting companion, the elliptical deformation can come from 
the static and periodic terms of the tidal potential as seen from 
the mantle frame of reference. In this case, periodic terms lead 
to tides, and the static term leads to the so-called static (tidal) 
bulge (the solid layer behaves as a fluid layer on the long term). 
Second, a so-called frozen bulge, resulting from previous states, 
may exist, as for instance in the Moon JGarrick-Bethell et al.l 
2006). In this case, the body is not in hydrostatic equilibrium. 
For our purpose, it is sufficient to distinguish between the re- 
sponse to tidal potential periodic terms, which we call the dy- 
namic tides, and a permanent (or very slowly changing) bulge, 
either owing to the tidal potential static terms or due to a frozen 
bulge, which we call a static bulge. The usual phase lag between 
the tide and the gravitational potential of the host body, which is 
due to internal dissipation, is not relevant for our purpose so is 
forgotten. Three dimensionless numbers are needed to describe 
the system: (i) the ellipticity p of the elliptical deformation, (ii) 
the Ekman number E - v/{Q. Rf), where v is the fluid kine- 
matic viscosity, R2 the outer radius of the rotating fluid, and Q 
its typical angular velocity before any instability, equal to the 
mean value of the (possibly varying) mantle spin rate Q.spi„{t), 
(iii) the differential rotation aQ between the fluid and the ellip- 
tical distortion, non-dimensionalized by the fluid rotation rate, 
aQ/Q. We distinguish two cases: a non-synchronized body, and 
a synchronized body. In the former, over one spin period, a mean 
differential rotation exists between the elliptical deformation and 
the fluid, whereas in the latter the mean rotation rates of the de- 
formation and of the fluid are equal. The different cases are de- 
scribed in the following and summarized in Table [1] 



D. Cebron et al.: Elliptical instability in terrestrial planets and moons 



2.1. Non-synchronized bodies 

For a non-synchronized body, we consider two cases depending 
on the origin of the elliptical shape. First, if the spin rate of the 
mantle Q.spi„ is constant, a TDEI can be excited by the tidal ellip- 
tical distortion due to dynamic tides, which rotate at the orbital 
velocity Q.„ i-h. This is th e standard configuration considered for 
instance bv lCraitd ( Il989l) . who showed that TDEI is indeed pos- 
sible except in a forbidden zone Q.spmlQ.orb £ [-1; 1/3], where 
no triadic resonance is possible. 

Second, if the elliptical shape comes from a static bulge, a 
non-zero mean differential rotation over one spin period implies 
that the fluid does not rotate at the same rate as the mantle, which 
corresponds to a spin-up or a spin-down process. This can occur, 
for instance, transiently after a large meteoritic impact, which 
is capable of fully desynchronizing the body (see for instance 
the considered scenario for explaining the Moon's magnetic field 
by Bars et al. 2011). In this case, a differential rotation exists 
between the fluid and the mantle with its static bulge, up to the 
typical spin-up/spin -down time neces sary for the fluid to recover 
the mantle velocity (iGreenspanll 1 9681) . i.e. up to 



t ■ - O ' F 

Ispin-up — ^'■spin ^ 



-1/2 



(1) 



Then, if the growth time of the TDEI is short enough compared 
to the spin-up time, one can expect a quasi-static evolution of 
the system, where the modification of the spin rate of the fluid 
is neglected during the growth of the instability: the former con- 
figuration is then transiently recovered. 

2.2. Synchronized bodies 

In the synchronized case, even if there is no mean differential 
rotation between the elliptical deformation and the fluid, oscilla- 
tions can nevertheless occur for different reasons. For the study 
of the elliptical instability, it is necessary to know the amplitude 
of these oscillations, which depends on their origins. We distin- 
guish below the forced librations caused by gravitational inter- 
actions with other celestial bodies, and free librations induced 
for instance by a meteoritic impact. 

2.2.1. Forced librations 

In forced librations, static bulge and dynamic tides have 
to b e considered simultaneously. To illustrate this, following 
iGol dreich & Mitchell (2010), we consider a simple toy model 
without any internal dissipation: a synchronously spinning satel- 
lite, with an elastic outer shell and a homogeneous fluid interior, 
moving on an elliptic orbit. The orbital velocity changes along 
the orbit, and writes at first order in the orbital eccentricity e as 



Q.orb = fi (1 + 2e COS Q.t), 



(2) 



where Q is the mean value of the mantle spin rate Q.spi„{t). 
Considering the influence of the orbital velocity variations on 
the satellite dynamics, we expect two limit cases: (i) if the rigid- 
ity of the elastic shell is zero or if the planet spin rate is low 
enough for the shape of the planet to have time to adapt to the 
gravitational constraints, the shell slides over the fluid and main- 
tains its equilibrium shape, with the long axis of the ellipsoidal 
figure pointing toward the companion body; (ii) if the rigidity 
of the elastic shell is strong enough or if the planet spin rate is 
rather high, the entire satellite rotates rigidly with a fixed shape. 
In the first case, only the elastic energy E^ias varies: the 
meridians of the shell are stretched and compressed due to the 



rotation, whereas the spin velocity of the satellite remains con- 
stant. This is the so-called optical libration. In the second case, 
only the gravitational energy Egrav varies and the spin velocity of 
the satellite changes, which corresponds to the so-called physi- 
cal libration. 

In both cases, the libration period remains small compared 
to the typical spin-up/spin-down time ([T), which means that the 
fluid does not follows the solid boundaries because it never has 
enough time to adapt to the periodic velocity fluctuations and 
continues to rotate at the constant synchronous rotation rate Q. 
This is the so-called no spin-up condition. In the first case, there- 
fore, a differential rotation exists between the fluid rotating at 
the constant rate Q. and the dynamic tides rotating at the oscil- 
lating orbital velocity, Af2/Q - 2ecosOf. An LDEI can thus 
be excited by this optical libration, as shown theoretically by 
iKersweU & MalkusI ( 119981) and iHerreman etall (2009). In the 
second, a differential rotation aQ/Q = e cos Qf exists between 
the fluid rotating at the constant rate Q and the static bulge sub- 
ject to physical librations of amplitude e, which depends on the 
internal structure of the satellite. The amplitude of the physi- 
cal librations e is always less than the 2e extreme value given 
by optical libration, because of different internal torques such as 
the gravitational torque and elastic strain torque (see for instance 
Van Hoolst et al. 2008, 2009). 

With the more general case of an a rbitrary torque applied to 
the shell, IGoldreich & Mitchelll (1201 Ol) estimated the ratio % = 

*^elasl ^grav ^Y 



% = 



32 n 1 -H V (1 -H kfY jldR^ 
~Y' 5 + v Ff GM^ 



(3) 



where v is the Poisson ratio, kf the fluid Love number, fi the sheU 
rigidity, R and d the mean radius and the t hickness of the shell, 
M the mass of the satellite. According to Goldreich & Mitchelll 
(1201 ol) . typical values give % ~ 10"^ for the subsurface ocean of 
Europa, and 51 ~ 0.1 for the subsurface ocean of Titan, whereas 
the silicate mantle of lo is expected to behave in the limit ^ » 1 . 
Because of the visco-elastic rheology of real bodies, the effec- 
tive response should be between the s e two extrem cases given 
by this model. iGoidreich & Mitchelll (1201 Ol) argue that the total 
increase in energy would be minimal, which leads us to consider 
that Europa and Tita n, for i nstance, behave l ike entirely fluid 



satellites. In cont rast. iKarate kin et al .l (l2008l), IVan Hoolst et al.l 
(I2008ll2009l) . and lBaland &Van Hoolst' ("2010^ consider that the 
arguments proposed by IGoldreich & Mitchell (2010) are unreal- 
istic on short timescales, hence that the rheology does not allow 
the bodies to reach their minimal energetic state. They assume 
that they behave rigidly, however, large libration may still be due 
to resonances with free libration modes, which may be reached 
for thin ice shells (Baland & Van Hoolst 2010). In addition to 
the displacements of the ice shell and mantle induced by gravi- 
tational interactions, relatively large longitudinal displacements 
may also be induced directl y in the fluid lay ers by the periodical 
part of the tidal potential (iTobie et al.ll2005l) . This would consti- 
tute a supplementary origin for a non-zero differential rotation 
between the fluid and the deformation. Either way, all these is- 
sues are still being debated and are clearly beyond the scope of 
the present paper All that is needed here is to know the ampli- 
tude of the relative motion between the elliptical distortion and 
the fluid. In the following, we consider the full range of con- 
figurations up to a maximum distortion amplitude given by the 
extreme value of optical libration. 
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2.2.2. Free librations 

After a meteoritic impact, for instance, so-called free libra- 
tions can occur on the typical resynchronization time (e.g. 
[Williams et al. 2001). Following the no spin-up condition ex- 
plained in section |2. 2. II the fluid keeps rotating at the orbital 
velocity (synchronized state), while the mantle librates around 
this mean value. The amplitude e of the free librations depends 
initially on the impact strength and decreases through time, and 
the libration frequency remains equal to a proper frequency of 
the body, given by Ufree - ^ V3(B - A)jC at first order in the 
orbital eccentricity, where {A,B,C) are the three principal mo- 
ments of inertia of the body (see for instance Lissauer 19851). 
Considering a static bulge, free librations can thus drive an LDEI 
from the differential rotation Afi/Q = ecos(w/„(.f)j providing 
that the growth time of the instability is shorter than the resyn- 
chronization time, as shown in section [375] 



Iso-V 

gravity anented by ^ = —V V 



"..(0 




Tidally deformed layer 

ellipticityp 



2.2.3. Zonal wind induced by pliysical librations 

Finally, in all scenarios involving phys ical lib r ation s, it has re- 
cently been determined analytically by Busse (2010), and con- 
firmed experimentally and numerically by Sa uret et al. (2010), 
that non-linearities in the Ekman layer driven by the librating 
rigid boundaries induce a differential rotation in the fluid of am- 
pHtude aQ/O = -0.154 (0coo/^)^, where 6 is the amplitude 
angle of the libration and ajg its frequency. TDEI can thus be ex- 
cited by this differential rotation with both static bulge and dy- 
namic tides. Nevertheless, the differential rotation generated by 
this process is always very small. We do not expect this mecha- 
nism to play an important role in a planetary context, since it is 
always dominated by LDEI, but it is worth here mentioning its 
existence since it may be relevant in certain astrophysical cases. 



2.3. Typical amplitudes of gravitational distortions 

The amplitude /3 of gravitational distortions, defined here as /3 - 



|flj - a^l/ia^ + a^), where oi and 02 are respectively the long and 
short axes of the outer boundary of the considered fluid layer, 
is generally unknown for celestial bodies. To study the elliptical 
instability for real cases, we need to estimate it, for instance by 
assuming an hydrostatic equilibrium shape. 

The equilibrium shape of a body of mass M and radius R, 
is an old pro blem that begins with the static bulge theory of 
iNewtonI (1686). This classical theory considers an incompress- 
ible no-spinning body at rest, deformed by a tidal field at leading 
order in R/D, which leads to a spheroidal shape and 



_3M2R^ 
^~2m"d3' 



(4) 



where M2 is the mass of the body responsible for the gravita- 
tional field and D the distance between the two bodies. This tide 
is sometimes referred to as the marine tide, where the gravita- 
tional potential of the tidal bulge is neglected. This approxima- 
tion always leads to a relevant but underestimated tidal deforma- 
tion. When possible, we use in the following a better estimate of 
j3 that takes the density distribution in the body and the gravita- 
tional potential of the tidal bulge into account: 



/3 



3 , M2R^ 

— "2 T 

2 M D^ 



(5) 



with the radial displacement Love number /12, directly linked to 
the potential Love number k2hy h2 - I + k2. A typical value is 



Fig. 2. Description of an internal liquid layer The tidal forces 
deform the core-mantle boundary (CMB) in an ellipse of axes 
fli and 02, which leads to an ellipticity /3 of the streamlines. In 
this layer, the fluid is rotating at the rate Q.. 



k2 = 3/2, obtained for an incompressible homogeneous body in 
hydrostatic equilibrium (e.g. Greff-Lefftz et al. 2005) . The tidal 
Love nu mbers can be calculated with the Clairaut-Radau theory 
(see e.g.|\MJiooLst et al."2008'). 

As shown by equation ([Sj, gravitational distortions vary with 
the interbody distance D. They can thus be divided into a compo- 
nent of constant amplitude, corresponding to the mean value of 
the gravitational distortion along the elliptic orbit, plus a smaller 
component with an amplitude oscillating between +3e times the 
constant one (e.g. Greenberg et al. 2003). In real cases that we 
consider in the following, this oscillating component can be ne- 
glected since it will always have a second-order influence on 
the elliptical instability compared to the constant (static or dy- 
namic) component of yS (but see appendix O. Nevertheless, in 
synchronized satellites, these so-called diurnal tides have impor- 
tant consequences for the internal state and the orbital evolution, 
since the changing shape of the bulge generates time-varying 
stresses, which generate heat by viscosity or friction. Besides, as 
shown bv'To bie et al.l (l2005l) . these diurnal tides should also in- 
duce longitudinal motions in the fluid layers, which constitute a 
supplementary origin of a non-zero differential rotation between 
the fluid and the deformation. As already mentioned above, in 
the context of the present study, this additional component is 
fully included in the amplitude of the considered physical libra- 
tion available to drive the elliptical instability. 

3. Generic formulas for the growth rate of the 
elliptical instability in a planetary context 

3.1. Model, equations, and dimensionless parameters 

We consider a telluric celestial body in the general framework 
sketched in figure [1] and we focus on a liquid layer described in 
figure |2] All dimensional parameters are listed in Table |2] The 
instantaneous spin rotation rate Q.spi„ - Q. (I + W(,)) may de- 
pend on time because of either free or forced physical librations. 
We focus on an internal fluid layer enclosed in an ellipsoidal 
shell, with an outer boundary of mean radius R2 at temperature 
62, and an inner boundary at temperature 61, with a mean radius 
Ri - r]R2. As already seen above, because of the no spin-up 
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Table 2. List of the dimensional variables used in this work. 



Table 3. List of relevant dimensionless parameters. 



Orbital rotation rate 


n<„-i(0 


Spin rotation rate 


n.p,„(f) 


Fluid rotation rate* 


Q. 


Mass of the deformed body 


M 


Mass of the attractor 


M2 


Inter-body distance 


D 


Free libration angular frequency 


UJfree 


Fluid layer mean external radius 


Ri 


Long/short axis in the equatorial plane 


Ol/flj 


Fluid layer mean internal radius 


Ri = r,R2 


Imposed external temperature 


O2 


Imposed internal temperature 


01 


Gravity at the external radius 


go 


Imposed magnetic field 


SoCxj 


Fluid density 


Po 


Fluid kinematic viscosity 


V 


Fluid thermal expansion 


1? 


Fluid thermal diffusivity 


V,h 


Fluid electrical conductivity 


O-e 


Fluid magnetic permeability 


t^ 



*mean value of the spin rate. 

condition, this fluid layer is initially rotating at the constant rate 
Q, equal to the mean value of Q.spin- This layer is considered to 
be homogeneous, with density po, kinematic viscosity v, ther- 
mal expansion 1?, thermal diffusivity v,/,, electrical conductivity 
(Te, and magnetic permeability yU. We focus here on the stability 
of the elliptical flow in the equatorial plane, but note that our 
local analysis remains valid in any plane orthogonal to the rota- 
tion axis. We choose R2 as the length scale and O ' as the time 
scale so that the mean basic spin of the body has an unit angular 
velocity along the rotation axis {O, tx,). The elliptical deforma- 
tion has a dimensionless angular velocity y^t) Cxj, with j^t) equal 
to Q.„rh{t)/0. when looking at dynamic tides and to Q.spi„(t)/Q. 
when looking at static bulges (see Table[T]). 

We consider the frame where the elliptical distortion is fixed, 
which is rotating at the angular velocity 7(1)^x3, with exj in the 
direction of the long axis oi and ex^ in the direction of the short 
axis fl2- The dimensionless equations of fluid motions are 



V -u = 0, 



(6) 



du dy(rt T 

— -H 2y ex3 X u H- -^ ex, x r h- u ■ Vu + V/9 = £ V^u + f , (7) 

where u is the fluid velocity; p the pressure (including the cen- 
trifugal term) non-dimensionalized by po Ri^ ^^', E = v/(Q/?2^) 
the Ekman number based on the external radius; and f = fs H- Tl 
the volumic force, including the buoyancy force fe and the 
magnetic Lorentz force t\^. The flow is rotating within an el- 



lipse Xj/flj + Xj 



/a: 



1, and we define the ellipticity as /3 



Using the dimensionless temperature - (6 - 02)/(^i - Oi), 
the temperature equation is 



de 

dt 



+ (u-v)e 



Pr 



(v^e - k). 



(8) 



where Pr — v/vti, is the thermal Prandtl number and K stands 
for a possible volumic heat source. Considering a gravity g = 
g(r,(ti)go^g, where e^ is a unit vector, (r, (f>) the cylindrical co- 
ordinates in the equatorial plane, and go the gravity at the ra- 
dius R2, the dimensionless buoyancy force to add in the Navier- 
Stokes equations using the Boussinesq approximation is f^ - 



Aspect ratio of the shell 


'7 


Ellipticity of the distortion 


y3 = (a2-fl2)/(fl2+4) 


Distortion rotation rate 


7(0 


Orbital eccentricity 


e 


Physical libration rate 


0)(t) 


Physical libration amplitude 


e 


Volumic heat source 


K 


Ekman number 


E = v/(Q R2^) 


Thermal Prandtl number 


Pr = v/v,i, 


Magnetic Prandtl number 


Pm = cr^ fi V 


Modified Rayleigh number 


Ra = &[et- e2]go/n^R: 


Magnetic Reynolds number 


Rm = 0-^11 QR2^ 


Elsasser number 


A = a-, Bl/(po fl) 



Ra 6 g(r,<i,) ^g, with the modified Rayleigh number Ra - § [61 — 
Oilgol^^Ri- In a planetary context, the temperature contrast to 
take into account only corresponds to the non-adiabatic compo- 
nent, which is the deviation from the thermodynamical equilib- 
rium state. 

We also take the possible presence of an uniform imposed 
magnetic field Bq along the rotation axis ex, into account, which 
is used as the magnetic field scale. The magnetohydrodynamic 
(MHD) equations then have to be solved simultaneously 



V B 







5B 1 , 

— -H (u ■ V) B = (B ■ V) u -H V^B 

ot Rm 



(9) 
(10) 



with the magnetic Reynolds number Rm - <Je ji QJi2 ■ The 
magnetic Lorentz force acting on the flow is given by Tl = 
(A/Rm) (VxB)xB, with the Elsasser number A = 0-^ Bq/(po ^)- 
All dimensionless parameters are Usted in Table |3] 



3.2. Base fields 

In the reference frame where the elliptical deformation is sta- 
tionary, the differential rotation of the fluid has an amplitude 
1 - 7(,). Besides this, the ellipticity induces an elongational flow 
-( 1 - 7(,)) p (x2 exi + xi exj), leading to the general elUptical base 
flow 



U = (1 - rw) [-(1 +y0)-^2exi + (1 -j3)xitx,] 



(11) 



This flow represents the laminar response of the fluid to the tidal 
distortion as an exact, non-linear solution of Navier-Stokes equa- 
tions for any finite viscosity, provided that (V x f) ■ es = 2 d,j{t). 
This means that a body volumic force noted f is necessary to 
have spin period fluctuations (i.e. free or forced physical libra- 
tions), as is obvious in a planetary context. Equation (fTTT i leads to 
elliptical streamlines of instantaneous ellipticity /3. We note that 
/3 is not the mathematical eccentricity of the streamlines, given 
by -y/2yS/(l +P). Also the velocity magnitude changes along a 
streamline, where the isovalues of the velocity are elliptical but 
with an ellipticity 2 /3. 

We further assume that a stationary temperature profile 
0(r, (f)) is imposed, which is, at order 1 in /?, the solution of 
the energy conservation equation (jS) with the base field U. We 
suppose that the modified Rayleigh number Ra is such that 
Ra = 0(J3). We also consider the presence of an imposed uniform 
magnetic field along the rotation axis, produced for instance by 
a companion body. We assume that the Lorentz force does not 
modify the base flow but only plays a role on the elliptical insta- 
bility. This implies that this force is 0(J3). In this context, we see 
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below that, regarding the elhptical instabiUty, equations for fluid 
motions on the order of in y6 are similar to those in the purely 
hydrodynamical case, where the elliptical instability is described 
as a resonance between two inertial waves. The magnetic and 
thermal fields only induce a correction in the fluid equations on 
the order of 1 in p, hence a correction of the growth rate of the 
instability, because of the stabilizing effect of the Lorentz and 
buoyancy forces. 

3.3. The WKB method: stability along a streamline 

Our local approach is base d on th e s hort-w av elength Lagrangian 



developed in iBavlvl (Il986h. ICraik & Criminalel 



then generahz ed in iFriedlander & VishiM (Il99ll) 



theory 

(Il98( 

and lLifschitz & Hameiril(ll99lh . This method has been success- 
fully applied to the elliptical instability by Le Dizes ( 200Q), 
then extended to take the energy equation a nd the buoyancy 
force into account (,Le Bars (fe Le D izesll2006l), or t h e indu ction 
equation and the Lorentz force [Herreman e t al.l (l2009l) . To 
summarize, the WKB (Wentzel-Kramers-Brillouin) method 
consists of looking for a perturbed solution of the equations 
of motion under the form of localized plane waves along the 
streamlines of the base flow. We thus look for a solution of the 
linearized non-dimensional system of equations (l6] [T] [8] |9] [TOl l 
in the form 



(12) 
(13) 
(14) 
(15) 

(16) 



U(x,,) = u+u'(0e■''»■^ 

P(.,o=P + p'(Oe"'*'^", 

B(,,o=Bo + b(f)e"'«■^ 

along the streamlines of the base flow described by 

At 

where k(,) is the time-dependent wave vector, x the position vec- 
tor, and where U (with its corresponding pressure field P), © and 
Bo = (0, 0, 1) are the dimensionless base fields defined in section 
13.21 Dropping the primes for simplicity, the linearized system of 
equations writes as 

k-u = (17) 

d,u -H i u (d,k ■ x) + i (U ■ k) u + (u ■ V) U + 2 y(,) Cxj x u 

= -i p k - yt^E u + — (i k X b) X Bo -H Tffl ? e„ (18) 
Rm 

A,e + lO (d,k • x) H- i (U ■ k) 6* -H (u • V) = -k^—O (19) 

Pr 

k b = 

d,b + ib(d,k-x) + i(U-k)b 



(20) 



(b ■ V) U + i (Bo ■ k) u - — b. 
Rm 



(21) 



Those equations can be decoupled in space and time to give 
an equation for the wave vector only: 



d,k ■ X -H U ■ k = 0. 



(22) 



The solution of the remaining equations for u, 6, b is then sought 
under the form of a Taylor expansion in p of all variables, as 
illustrated in appendices lAl and iB] 

This approach is used in the following sections to calculate 
the growth rate of the instability in the two generic cases: the 
TDEI, which appears in the case of non-synchronized bodies, 
and the LDEI, which appears in the case of synchronized bodies. 



3.4. Non-synchronized bodies: inviscid growth rate of the 
TDEI 

In this section, we consider the effects of dynamic tides of am- 
plitude j6 on the liquid core of a Mercury-like planet orbiting 
close to its star with (i) constant but dilTerent orbital and sidereal 
rota tion periods, (ii) an imposed thermal stratification (see e.g. 
Ma nglik et al.ll2010h . and (iii) an externally imposed magnetic 
field (e.g. the Sun magnetic field). The same analysis applies to 
the stratified zone of a star (the so-called radiative zone) tidally 
deformed by a companion body, taking the magnetic field gener- 
ated by dynamo in its convective zone into account. The present 
configuration corresponds to the standard case of the elliptical 
instability as already known, but completed by the complexi- 
ties present in real astrophysical cases. These additive effects 
have already been studied separately, even if they are simulta- 
neously present in real systems. The effect of the angular veloc- 
ity of the tidal bulge has been studied in Mivazaki & Fukumota 
( 199|). lLe DizesI (l2000h . lLe Bars et all (l2007h and lLe Bars et all 
(.2010|) and the presen c e of a thermal field has been studied 
in iLe Bars & Le DizisI (l2006h . who study the linear competi- 
tion between the growths of the TDEI and the convection, as 
well as in Cebron et al. (20103) and|Lavorel & Le Bars (2010), 
who study the growth of the TDEI over established convec- 
tive flows. The prese nce o f an inner solid core has been stud- 
ie d in lUa caze et alj (l2005h an d of an external ma gnetic field 
in 'Kerswell ( 1994), 'Kerswel]| ( l2002h . iLacaze etalj (12006^ and 
Herrem an et al. (2009) . We extend these works by including all 
of these features in a single formula. 

In the non-synchronized case, which is considered in this 
section, the base flow (fTTI) reduces to 



U = (1 - Qorb/^spin) [-(1 +P)x2exi + (1 - /S)xie^J. 



(23) 



The WKB analysis is then tractable (see appendix [A|, taking 
thermal and magnetic effects into account in the limit where 
buoyancy and Lorentz forces are on the order of j8. The insta- 
bility does not exist in the range Q.spi„/D.„rb e [-1; 1/3], which 
is called the forbidden band. It corresponds to the absence of 
resonance between the elliptical forci ng and the inertial waves 
of the rotating flow (see lLe Bars et al.li2007.. for a complete dis- 
cussion). In the present limit, the presence of the thermal and 
magnetic fields does not affect the forbidden band. Neglecting 
the thermal diffusion, the inviscid growth rate of the TDEI with 
the presence of thermal and magnetic fields is 



-y/(2QG + 3)4 /32 - 4 (Ra r d,.®f 
16 11 +OG|3 



A 



4|1 +QC|3 



(24) 



with Q.'-' - Qorb/i^spin - ^orb), '" the radius and drd the dimen- 
sionless temperature base-field radial gradient on the considered 
streamline. This expression allows to recover the different cases 
already obtained in the literature. For inst ance, the purely hy - 
drodynamic gro wth r ate given by|M ivazaki & Fukumotol(ll992l) . 
iLe DizesI (l20()0l) . and lLe Bars et al.i IQQI 0) is recovered for (Ra = 
0, A = 0). For a fixed elhptical deformation, we recover the clas- 
sical inviscid value o-i„t,//3 - 9/16. Finally, in the absence of a 
thermal field and with a stationary bulge (Q.nrh - 0, Ra = 0), the 
magnetic case given in lHerreman et al.l (l2009l) is also recovered. 
Formula (l24l l is fully generic and clearly illustrates the stabiliz- 
ing influence of Joule dissipation and of a local stratification in 
the range of validity of this stability analysis (see section lXTI and 
appendixlDt. 
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3.5. Synchronized bodies: inviscid growth rate of the LDEI 

In this section, we consider the Uquid elUpsoidal core of a syn- 
chronized moon Hke lo or of an extrasolar telluric planet or- 
biting close to its massive attractor with (i) an orbital period 
equal to the sidereal rotation period, but with small instanta- 
neous fluctuations of the differential rotation between the el- 
liptical deformation and the fluid, whatever their origin (opti- 
cal or physical, forced or free libration, longitudinal flows in- 
duced by tides); (ii) an imposed magnetic field (e.g. Jupiter's 
magnetic field for lo); and (iii) a local thermal gradient. The di- 
mensionless instantaneous differential rotation between the fluid 
and the elliptic deformation oscillates with a libration ampli- 
tude e (equal to 2e for optical librations) and a libration fre- 
quency ojo (equal to 1 for forced librations). Considering the 
particula r case of a fluctuation due to th e orbital ellipticity (i.e . 
oJo = n. lKerswell & MalkusI (Il998h and lHerreman et alj (l2009l) 
have shown that these oscillations can lead to LDEI. We ex- 
tend these previous studies to the more general case of small 
fluctuations of arbitrary periods, taking buoyancy into account. 
Neglecting the thermal diffusion, the inviscid growth rate of the 
LDEI is then (see appendixlBJ) 



16 + ajl 
64 



U._^(K„.,,8)=_^A 



(25) 



at first order in e/3, taking the effects of thermal and magnetic 
fields into account in the limit where buoyancy and Lorentz 
forces are 0(J3). The forbidden band is given by \a>o\ > 4. As 
before, the generic formula (l25t clearly illustrates the stabilizing 
influence of Joule dissipation and of a local stratification. 

The case oj,, — > corresponds to the limit toward the fully 
synchronized state. In the case of the TDEI, this limit case is 
obtained with Qp — > oo, which gives the inviscid growth rate 
Cinv - y6/(4Q'^) for large wavenumbers. Both expressions for the 
growth rate are thus consistent in the limit of synchronized state: 
Q.orb/Q.s i,in = i - e i.e. | 1 -h Q.^ \ ~ |Q*^| = 1/e. The expression 
given in iHerreman et al.l (l2009l) is also exactly recovered when 
(L),, - 1 and A = 0. There is a slight error on t h e num erator of 
the magnetic damping term in IHerreman et al.l (1200 9*): in their 
considered case, aio = 1, the magnetic damping is erroneously 
-A/16, instead of -3A/ 16. 

3.6. Viscous dissipation 

The previous sections present the calculation of the growth rate 
of the elliptic instability in an inviscid fluid with Joule dissipa- 
tion and buoyancy stabilization. Calculation of the threshold of 
the instability requires correctly estimating all dissipative terms. 
In the case of no-slip boundaries, dissipation occurs mainly in 
the viscous boundary layers of thickness E^^~. This implies a 
damping term that should alter the growth rate: 



o-^o-iin.-aE""- f(7]). 



(26) 



where a is a constant between 1 and 10, equal to a = 2.62 and 
firf) = (1 + r]'^)/(l - rf') for the spinover mode of th e TDEI (see 
e.g. Kudlick 1966; Hoflerbach & Kerswell 1995; Lacaze et all 
120051) . 

In addition to decreasing the growth rate, viscous dissipation 
is also primordial for quantifying the orbital evolution and rota- 
tional history of a binary system during its synchronization. A 
model has been proposed in Le Bars et al. (2010) for Q.orb - 0, 
which allows the authors to estimate the viscous power dissi- 
pated by TDEI. Our purpose here is to generalize this model to 



all cases studied a bove. Far from threshold, the model proposed 
by 'Le Ba rs et al.l (1201 Ol) considers that the TDEI simply corre- 
sponds to a differential rotation between the boundary and the 
bulk. According to this model, the power dissipated by the sys- 
tem is 



P= -2M/?? AO^Q£''^^ 



(27) 



assuming that in the small Ekman numbers limit reached in as- 
trophysical cases, the amplitude of th e instability is comm ensu- 
rate with the differential rotation aQ dCebron et al.ll2010al) . 

The tidal quality factor Q is widely used in systems evolution 
calculations. By analog y with the theo ry of harmonic oscillators, 
Q is defined by (e.g. Greenberg 2009, for a recent discussion 
on Q) the ratio between the maximum potential gravitational 
energy stored in the tidal distortion over the energy dissipated 
in one work cycle. The dissipated power associated to the flow 
driven by the elliptical instability does not have the periodicity 
of the forcing, so a quality factor cannot be rigorously defined 
in the same way. However, we can define a closely related ratio 
Q* , comparing the power typically dissipated in the fluid layer 
over one revolution with the potential energy stored in the bulge, 
which is on the order of £■(> ~ 47rpo ^p s^ P S,, with s the dimen- 
sional height of the tides (e.g. lBenesil 19901) . Since /? ~ SIR2, we 
obtain from equation ( l27l i 



Q* 



80/^^ 



/?AQ2£l/2' 



(28) 



The ratio Q* gives a dimensionless measure of the strength of 
the dissipation in the fluid. 

3. 7. Validity of the approach 

The previous analysis is valid when the elliptic instability 
comes from a resonance of pure hydrodynamic inertial waves. 
Therefore, any previously derived expressions are limited to the 
case where buoyancy and Lorentz forces are Oifi). According to 
( ITTI i. this means that 



k ~/3 and Ra ~ j3, 

Rm 



(29) 



where k is the dimensionless norm of the wavevector of the ex- 
cited mode. For a typical planetary core, these conditions can be 
rewritten 



Bo ~ 0.1 



'/3/k 



RiE 



in yuT and 



non-adia 
^ adia 



10 



-3 



/3 



E'Rtsl' 



(30) 



where F„„„-adia and Fadia are the non-adiabatic and adiabatic 
components of the core heat flux, respectively. The condition 
on the magnetic field is easily verified for planetary cores over a 
wide range of wave vector k. The condition on the non-adiabatic 
heat flux is more problematic to quantify: in most planets, the 
adiabatic profile is supposed to be sufficient to transport core 
heat flux, and the non-adiabatic component is estimated to be 
very small, but it is not known precisely. One should notice that 
the condition ( l30l l is very restrictive, and special attention should 
be paid in each given configuration. For instance, in the case of 
Europa, Eq. (l30l l implies F „on-adia I F adia ~ 0.1%, which seems 
reasonable; but in the case of lo, it implies F„„„-adialFadia ~ 
1%, which is only marginally verified since the estimated non- 
adia batic heat flux is about one fifth of the adiabatic component 
(e.g. lKersweU & Malkusl[T998h . 
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Table 4. Physical and orbital characteristics used for the stability 
calculations and results. 



ir 



MCxlO^-^kg) 
R(km) 

-* spin \W/ 
Torb (d) 

Tidal amplitude (m) 

Ri/R 

V 

p (measured) (xlO^) 
j8 (hydrostatic) (xlO^) 
B.,/ (nT) 
Bo' (nT) 
A 



Mercury 

0.330 

2440 

58.6 

87.97 

0.925° 

0.8" 



21 

7.6 

6.8 

250'^' 

488 

6.4 • lO^*^ 

-1.5- 10-"^ 



Venus 

4.87 

6051 

-243 

224.7 

1.8 

0.17 



316 

? 

1.1 

30 
6100 
0.004 
-5.43 



Earth 

5.98 

6378 

0.997 

27.32 

0.6 

0.55 

0.35 

0.11 

1.9 

0.8 

3-10* 

1.8-10' 

0.015 

-7.7 



Early Earth* 
5.98 
6378 
0.418 
9.67* 
3.4 
0.55 



0.047 

11 

6.7 







0.003 



Following iHerreman et al.l ( 120091) . we take as 
cr^ = 4 • 10' S.m"', po = 12 000 kg.m"' and v 
consistent with a Fe/Fe-S composition. 



typical 
= 10-'' 



values 

^2 c-l 



" IVanHoolstetal J ( 120071) 

* Considering an Early Moon two times closer than today. 

^ Equatorial surface fiel d 

'^ I Anderson etal.l (1201 Oh 

*■ Considering a variation in r^^ from the core to the planetary surface 

(r being the spherical radius). 



Now supposing that the buoyancy or the Lorentz force is 
on the order of in yS, we can wonder if the elliptic insta- 
bility still exists. In this case, inertial waves are replaced by 
gravito-inertial or magneto-inertial waves, and the elliptical in- 
stability arises as a resonance between those modified waves. 
Resonances of magne to-inertial waves has been studied, for 
instance, in |K erswell (1994]), iLebovitz & Zweibell (12004) and 
iMizerski & Baler (.2009.) in the case of an imposed uniform mag- 
netic field along the spin axis. This is the so-called magneto- 
elliptic instability. Resonance of gravito-inertial wav es has 
been studied in Le Bars & Le Dizes (2006) and Guimba rd et al.l 
(l2010l) . who concluded that a stratified field can either be stabi- 
lizing or destabilizing depending on the shape of gravitational 
iso-potentials and isotherms: the so-called gravito-elliptic insta- 
bility. This point is further clarified in appendix iP] which shows 
the high sensitivity of the elliptical instability to the specifici- 
ties of the thermal and gravity fields. Planets with buoyancy or 
Lorentz force on the order of in y6 should be the subject of spe- 
cific studies, which is beyond the scope of the generic results 
presented here. 



4. Application to solar/extrasolar systems 

Using previous results, we are now in a position to calculate the 
threshold of the elliptical instability for telluric bodies of differ- 
ent systems. We consider that a body is stable or unstable when 
the mean value of the growth rate over an orbit is either positive 
or negative. 



4.1. Non-synchronized system: tide-driven elliptical instability 
(TDEI) 

We consider first the TDEI in liquid cores of telluric bodies in 
the solar system. A rough criteria given by the equations (l24T i 
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Fig. 3. TDEI stability diagram for celestial bodies of the so- 
lar system. Considering a surfacic viscous damping term of the 
growth rate cts - -a f{rj) E^^^ (see section [J!6T l. the zone above 
the black line, defined by ( l33l l. is the unstable zone, whereas 
y, defined by (ISTT i. is calculated with a - \ and a - \Q for 
each planet (for the actual Earth and Venus, the difference is very 
small). The yellow zone is the so-called 'forbidden zone', given 
byQ,p;„/no,&e[-l;l/3]. 



and ( l26l l leads to a threshold pj y[E ~ 0(1), as already men- 
tioned. In the solar system, this leads to focus only on Mercury, 
Venus, and on the Earth-Moon system during its evolution. Most 
tidal evolution models predict that the Moon rapidly retreats to 
25 - 35 Earth radii in less than about 100 Ma (IWebblll982l: 
iRoss & SchuberJIl 9891: IWilliamsll2000l [2004 . and we thus con- 
sider two limit cases: the actual Earth-Moon system and an early 
Earth with an early Moon at 30 Earth radii, i.e. two times closer 
than today. The tabulated values found in the literature for these 
planets are given in Table |4] The Ekman number E is calculated 
with a molecular kinematic viscosity v = 10"^ m^.s"', consistent 
with Fe/Fe-S composition of a liquid outer core. 

We first neglect thermal effects. To represent all bodies of 
Table |4]on the same stability diagram, we define the quantity 



Y^pLfiri)^- 



A 



4|1 +QC|3 



(31) 



The quantity Y includes the specific dependence of the growth 
rate on the spin/orbit angular velocity ratio, aspect ratio 77 of the 
inner core and the magnetic field. The stability criterium 



(T - 



(20"^ 
16 



l^ + 3f ,- A 



derived from (l24l and (l26t is then equivalent to 



Y> 



16 11 -hQ' 



G|3 



(3 + 2nc)2 



(32) 



(33) 



Figure[3]represents the stability results for the TDEI in the liquid 
cores of the non-synchronized planets of the solar system con- 
sidered in Table |4] The case of each planet is discussed in the 
following. 

An important result, akeady noticed in'Cebr on et al.l(l2010ah . 
is that the Early Earth, with a Moon two times closer than to- 
day and in the absence of an external magnetic field, is unsta- 
ble with a good level of confidence. The dissipated power due 
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Table 5. Physical and orbital characteristics used for the four 
Galilean moons and Titan. 



20 30 40 50 

Earth - Moor distance ir Earth radii 

Fig. 4. Evolution of the typical growth time T growth - 1 1 cr of 
the instability in the Early Earth core as a function of the Earth- 
Moon separation and of the length of day on the Early Earth 
Tspin (o; - 2.62). This diagram assumes no external magnetic 
field (Bo = 0) and no thermal field {Ra = 0). 



to the instability was around 5 ■ lO'^ W, which corresponds to 
Q* ~ 0.003. This estimation seems huge in comparison to the 
presen t dissipation by tidal fr iction (~ 3.75 ■ lO'^ W accord- 
ing to ' Munk & Wunschlll998h . However, one must notice that 
the estimations given here are based on a model where the rota- 
tion rate is implicitly expected to be constant, corresponding to a 
quasi-static approximation of the orbital evolution of the system. 
This approximation obviously breaks down for a high dissipa- 
tion rate, and the above result should be interpreted as proof of 
a rapid orbital evolution of the Earth-Moon system. The Moon 
orbit inclination is not taken into account in the stability analy- 
sis considered in this work. A similar analysis would be diflicult 
because the forced base flow in such a configuration is not ana- 
lytically known. However, it has been demonstrated numerically 
in|Cebron et al. (2010a) that this inclination only has minor con- 
sequences on the process of elliptical instability, which may be 
easily taken into account in considering an effective elliptical 
distortion in the equatorial plane. In the present case, this orbital 
inclination would slightly decrease our growth rates but does not 
change the orders of magnitude our conclusions either 

Figure |4] shows the stability of the Early Earth in more 
detail, for different values for the length of day and Earth- 
Moon distances. In the absence of meteoroid impacts, the 
angular momentum conservation links these two quantities, 
which cannot vary independently. However, at this epoch, vi- 
olent meteoroid impacts have probably modifi ed the angular 
momentum of the early Ear th-Moon system (iMeloshI 119751: 
IWieczorek & Le Feu vrel 120091) . so we keep these two parame- 
ters independent, which allows us to cope with uncertainties. 

The case of the actual Earth is more subtle because if we 
consider that the actual magnetic field is provided by thermo- 
solutal convective motions in the core, it has to be considered as 
an imposed field for the dynamics of the elliptical instability. In 
this case, the destabilizing term in the growth rate (l32t is about 
10"'', whereas the magnetic damping term is around 0.004. Then 
the TDEI cannot grow, regardeless of the Ekman number. On the 
contrary, if we consider that the actual magnetic field is provided 
by the flow driven by the TDEI, the threshold has to be calcu- 
lated with A = 0, and the actual Earth is slightly unstable, with 
a growth time of around 14,000 years. The same result was sug- 
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the elliptical instability frorn grow ing, as shown in lCebronet al.l 
: lLavorel & Le Bars! (*2010). Since this temperature con- 



gested bv lAldridge et al.l (Il997h . but neglecting the influence of 
magnetic field and global rotation. 

We now consider the influence of thermal effects. 
Considering the actual heat flux of Earth, 'Chris tensen & AubertI 
(2006) estimate the corresponding super-adiabatic temperature 
contrast to be about 1 mK, leading to a vigorous convection 
in the Hquid core. Such a vigorous convection do es not prevent 
the ellip ti 
(I2010cl) :[ 

trast is uncertain, we can consider as an extreme case for stabi- 
lization of the TDEI, a subadiabatic gradient on the same order 
of magnitude. For the actual Earth, such a stratification leads to 
Ra « -1.4 ■ 10"^, and the dependence Ra oc D^ E^ gives the 
value Ra X ~2 ■ 10"^ for the Early Earth (see lSumita & Yoshidal 
. 2003. for a discussion of this possible stable density stratifi- 
cation in the whole early Earth outer core and its disruption). 
Considering Ra ^ -1 ■ 10"^ as an upper bound, formula (l24l l 
is valid (Ra ~ 0(J3)) and shows that the thermal stratification 
reduces the growth rate by 2%. This confirms that the role of 
the temperature can be neglected in the limits considered in this 
work. 

For the last two bodies, Venus is in the forbidden band, 
which means that no matter what the tidal deformation or the 
Ekman number is, the TDEI cannot grow. Mercury is slightly 
below the threshold of the instability and is thus probably sta- 
ble today. Mars is clearly stab le nowadays (0 «: V e), but 
lArkani-Hamed et al.l (l2008h and lArkani-Ha med (2009) suggest 
that past gravitational interactions with asteroids could have ex- 
cited a TDEI in the martian core during their fall towards the 
planet. 

4.2. Synchronized body: libration-driven elliptical instability 
(LDEI) 

4.2.1. Galilean moons and Titan 

The presence of LDEI in l o has been fir s t sug gested by 
iKersweU & MalkusI (Il998h . In lHerremanetalJ (12009) . this sug- 
gestion is reexamined and the magnetic field induced by this 
possible instability quantified. In the following, this calculation 
is re-evaluated and extended to the four Galilean moons (lo, 
Europa, Ganymede, and Callisto), considering the presence of 
the external magnetic field of Jupiter Titan is also considered. 
All necessary data are given in Tables|5]and|6] 

As described in section 13.51 we consider an instantaneous 
differential rotation ecos(wof) for these synchronized bodies. 
Focusing on the forced librations due to the orbital eccentricity, 
the libration frequency is w,, = 1 . The amplitude of the libration 
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Table 6. Stability results in the extreme case of optical librations 
(e = 2e and/? from this Table and Table |5]l. 



lo Europa Ganymede Callisto Titan 



RilR 

1 

£(xlO"') 

A(xlO^) 
cr Cyr"') 



Crust (km) 
Depth (km) 
RilR 

j8(xl0'') 
A(xlO''') 



core 
0.52"* 



2.7 

60" 

42 
0.016 



core 

0.38* 



14 

9.1'' 

4.1 

0.0025 

ocean 
10* 
100« 
0.99 
0.94 
2.0 
9.T' 
21 

0.0016 



core 

0.27'' 



20 

3.T 

0.7 

-3 ■ 10-* 

ocean 
100* 
150'' 
0.96 
0.94 
1.5 
3.7" 
3.5 
-6 ■ 10-* 



ocean 

150' 

150' 

0.94 

0.93 

4.5 
0.72'- 

0.9 
-4 • lO^* 



ocean 
70^ 
200^' 
0.97 
0.92 
3.5 
1.2/ 

-lO^* 



For the liquid cores, we take 0-^=4- 10^ S.m"', po = 8 000 kg.m"^ 
and V = 10"' m^.s"' as typical values, consistent with a Fe/Fe-S 
composition. For the subsurface oceans, we take cTe = 0.25 S.m"' 
( lHaiid&Chvb3 120071) . po = 1000 kg.m"^ and v = lO"*^ m^.s"' as 
typical values. 



Kerswell & Malku: 



Hussmann & Spoh 
Bland e t al. (2008) 



considering the static tidal bulge 



' with k2 « 0.3 dWahr et al.ll200g : lBaland & Van HoolslllMol) 

■■ Eq. (O with k2 a 03 

f k2^ 1 ( Goldreich & Mitchellll20Tol) 
^I Wahretal .1(12006) 
'' bland et all J2009 '^ 



'iKuskov & KronrodI Ooof) 
^ ISohl et al.. (2003) 



e is given by e = 2e for optical librations. For physical libra- 
tions, obtained for 51 » 1, e has to be measured or estimated, 
but are less than the extreme value 2e (see the data in lNoir et al.l 
[2009 ). The theoretical analysis is the same, and the use of for- 
mula (l25T l with Ra = Q (thermal field negligible) and ajg = 1 
gives the LDEI threshold. To obtain a unique stability diagram 
for all bodies in Table |5] we define the quantity 



Yi 



4 

£B A 

^ 17 



[a(l-7,)/(„)f 



(34) 



and use the Ekman number based on the thickness Ek - E/(l - 
T])^. The threshold for LDEI given by formulas (l25t and (l26t 

cr^^eP-a(l-7^)f^r,■>^|F,-^A>0 (35) 

is then equivalent to 

64 r- 
Y2>—Wk- (36) 

This allows us to plot the stability diagram shown in figure 
|5] in the extreme case of optical libration (e = 2e) for a 
quasi-equilibrium hydrostatic bulge calculated with formula (|5]), 
corresponding to the optimal case for LDEI (i.e. the maxi- 
mum possible libration amplitude and the maximum possible 
elliptical deformation). In the following, we discuss the stabil- 
ity versus the LDEI of the Galilean moons. Titan, and three 
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Fig. 5. LDEI stability diagram for synchronized celestial bod- 
ies in the optimal case for the instability, i.e. an optical libration 
(e = 2e) for a quasi-equilibrium tide. All values are given in 
Table |5] and the horizontal axis represents the Ekman number 
Ek based on the thickness of the fluid layer The label KM98 
for lo calls that the po int is placed with the values used by 
iKerswell & Malkiii (11998). The zone below the black line, de- 
fined by the viscous surfacic damping coefficient a - 2.62, is 
the stable zone, whereas the black dashed lines represent the ex- 
tremum values a —\ and a = 10. 



Super-Earths. Because the tidal bulge and the libration ampli- 
tudes of these bodies are no t yet known (see the discussion of 
iGoldreic h & Mitchelll |201()I) . we present the results in figures 
(|6]l-([8]) on diagrams in the (jS, e) plane, taking the full range of 
variability of /? and e into account. Therefore, the upper right 
hand corner will correspond to the optimal case for the LDEI: 
the libration of a purely deformable body i.e optical libration 
with e - 2e and an hydrostatic bulge. The lower left hand cor- 
ner coiTesponds to the libration of a rigid body (physical libra- 
tions), associated with the weak diurnal tides. In the same way, 
the lower right hand corner corresponds also to physical libra- 
tions, but with a hydrostatic bulge. Finally, the upper left hand 
corner corresponds to the libration of a purely deformable body 
(optical libration) associated to the small diurnal tides ampli- 
tude. The relevant physical configurations for each body depend 
on their compositions so is specifically discussed in the follow- 
ing for each of them. 

First, we consider lo wi th th e values used in the s tudies 
of Kers well & MalkusI ( Il99g ) and He rremanet al.1 (|2009), i.e. a 
static bulge of ellipticity yS = 0.006 and a libration amplitude 
assumed to be e = 2e (see Table |5]|. As already found by these 
authors, lo is unstable with a good level of confidence (typical 
growth time of 63 years). However, this is an optimal unrealistic 
case, because the ellipticity used is due to a static bulge (^ » 1), 
and the libration amplitude is taken as equal to 2e, as in the de- 
formable case (iR « 1). Due to its silicate mantle, the core of 
lo is expected to be in the limit 51 » 1, and consequently the 
ellipticity to consider is indeed /? - 0.006 but the libration am- 
plitude is instead e ^ 1.3 ■ 10"^ (see Table |5]l, which is 63 times 
smaller than 2e. With these more realistic values and condidering 
the presence of Jupiter magnetic field, lo is expected to be sta- 
ble, unlike what was expected. However, to obtain a better view 
of the stability in lo, figure |6] gives the typical growth time of 
the instability for different ellipticities and libration amplitudes, 
ranging between the limit cases 51 » 1 and 51 <K 1, with e = 2e 
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mo' T^j^j^tepaes(colorbar:10") 
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Fig. 6. Evolution of the typical growth time T growth - 1 /cr of the 
instability in lo with the tidal bulge ellipticity yS and the libration 
amplitude e {a - 2.62). The up per right corner cor responds to 
the optimal cas e considered in Kerswell & MaUciisI (1199 8) and 
iHerreman et al.l (12009?) . the upper left corner corresponds to op- 
tical Ubrations with a small tidal amplitude corresponding to di- 
urnal tides OS = 7 ■ 10"^ e = 2e = 0.0082) and the lower 
right corner corresponds to the region of physical Ubrations of 
the static bulge (JS = 0.006, e = 1.3 ■ 10""*). The white zone 
corresponds to the stable zone where the instability cannot grow 
because of dissipative effects (at the boundary with the colored 
zone, the growth time is infinite). The colorbar range is chosen 
so that color variations are visible. 



and the diurnal tidal ellipticity 3ey6 ~ 7 ■ 10"^, corresponding to a 
diurnal tide amplitude of 130 m. Future accurate measurements 
of the tidal amplitude at the core-mantle boundary and of the 
libration amplitude should confirm our prediction. 

In the optimal case for instability, figure |5] shows that the 
liquid core o f Europa is unstable, as already suggested by 
iKerswell & M alkus ( 1998), even when taking the Joule dissipa- 
tion due to the presence of Jupiter magnetic field into account. 
The typical growth time of the instability to be around 400 years, 
and the associated dissipation to be on the order of P ~ 10"^ W. 
This corresponds to Q* ~ 10^, and is two orders of magnitude 
below the conservat ive estimation of 3 ■ 10'^ W for the tidal heat- 
ing rate on Europa (lO'Brien et al.ll2002b . In reality, the silicate 
mantle of Europa should behave more rigidly. Since the libration 
amplitude of the mantle and the amplitude of the tidal distortion 
at the core-mantle boundary are not known yet, all intermedi- 
ate behaviors are explored in figure |7^. We conclude that high 
libration amplitude and/or rather large elliptical distortion are 
needed for Europa's core to be unstable; nevertheless, the pa- 
rameter range for instability is rather wide in figure ^ and seems 
to be reachable: we thus expect the core of Europa to be unsta- 
ble. Results for Europa subsurface ocean are shown in figure 05. 
The elastic behavior of the icy crust above the subsurface ocean 
is expected to behave in the deformable limit {% <s 1), which 
correspond to parameters close to the upper right hand corner. 
As also shown in figure|5] Europa's ocean is therefore unstable. 

Concerning the two last Galilean moons, the core of 
Ganymede and the subsurface oceans of Callisto and Ganymede 
are found to be stable in the optimal case (fig. |5]l. An LDEI is 
improbable today. Figure |5]also shows that the subsurface ocean 
of Titan is probably stable, because even in the optimal case for 
the LDEI, it remains in the vicinity of the threshold. 
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Fig. 7. Same as figure |6]but for Europa, considering (a) its pos- 
sible liquid core and (b) a 100 km depth subsurface ocean. Both 
the core and the subsurface ocean of Europa are expected to be- 
have in the deformable limit, locating their states in the upper 
right comer of the diagram, so both are found to be unstable. 



4.2.2. Super-Earths 

The recent discovery of extrasolar telluric planets gives typical 
examples of synchronized planets in close orbit around their host 
stars. This particular astrophysical configuration should lead to 
a vigorous LDEI in their possible liquid cores. In this sec- 
tion, we consider three Super-Earths, expected to be telluric: 55 
CnC e, CoRoT-7b, and GJ 1214b, at D = 0.0156 A.U, D = 
0.0172 A.U, and D = 0.0143 A. U from their host star respec- 
tively. T he data u sed here are from lWinn et a l.'('201 1*) for 55 CnC 
e, Valen cia et al.l (|2010) for CoroT-7b, and Charbo nneau et al] 
(2009) for GJ 1214b. They are given in Table Q For CoRoT-7b, 
the work of iLeger et al.l (l201ll) predicts a core composed of liq- 
uid metal, representing 11% of the total planetary volume, as for 
the actual Earth. Th e presence of a liqui d core in Super-Earths 
is still not clear (e.g. lWaener et al.ll20I iL for CoroT-7b), but we 
can reasonably assume that a planetary liquid core could occupy 
about one third of the planet radius, which corresponds to 4 % of 
the total planetary volume. Because of the proximity of the par- 
ent star, these extrasolar planets are expected to be synchronized. 
The actual orbital eccentricities of CoRoT-7b and GJ 1214b are 
not known. If they are fully circularized and synchronized, no 
elliptical instability can grow. In contrast, if a small libration 
exists, previous stability formula can be used. We assume here 
an orbital eccentricity of e = 0.001, well beyond the detection 
limit. Figure |5] shows that in the optimal case these three Super- 
Earths cores are clearly unstable with a good level of confidence, 
which means that the LDEI is probably present in their liquid 
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Table 7. Physical and orbital characteristics used for the stability 
calculations in exoplanets. 

CoRoT-7b GJ 1214b 55 CnC e 

M (in Earth's mass) 4.8 6.55 8.57 

/? (in Earth's radius) 1.68 2.678 1.63 

Torb(d) 0.854 1.58 0.7365 

Ri/R 1/3° 1/3* 1/3'' 

e O.OOr O.OOK 0.057'' 

7/ 

£(xl0"') 9.4 6.8 8.6 

j8(xl03) 7 6 5 

(T(yr-') 0.01 0.005 0.45 

We take as typical value v = 10"* m^ s^\ consistent with a Fe/Fe-S 
composition. We use the formula ^ to estimate the tidal bulge 
ellipticity. 



a T ^ in years (colorbsr: 10^ 



" Coherent with lLeger et a l. (2011). 

* Coherent with the values from iNettelmann et al.l (1201 ih . 

■^ Assumed hypothe sis. 

'' IWinnetalJ ( l20TIh . 



layers. One can see that this result is not very sensitive to the hy- 
pothesis on the size of the considered liquid core: for instance, 
with the values of Table|7] CoroT-7b becomes stable for a liquid 
core aspect ratio R2/R below 1 % (for the optimal case i.e. opti- 
cal librations and equilibrium tides). In these estimates, we use 
the hydrostatic tidal deformation, which underestimates the real 
tidal deformation. Also, higher orbital eccentricity would lead 
to more unstable configurations. Figure |8] shows the influence 
of these uncertainties, as well as the effect of smaller libration 
amplitudes and tidal deformations. The proximity of the known 
Super-Earths with their host stars leads to strong tidal deforma- 
tions, and the LDEI is then able to grow from very small libra- 
tion amplitudes, as shown in figure[8](the true tidal deformation, 
larger than the hydrostatic value, leads to an LDEI for a large 
range of libration amplitudes). Same conclusions are obtained 
for Kepler-lOb (Ba talhaet al. 20 1 1), assuming the same hypoth- 
esis. To conclude, the presence of the LDEI in their liquid cores 
is very probable, whatever the uncertainty ranges. 



5. Conclusion and discussion 

In conclusion, we have investigated theoretically the elliptical in- 
stability in telluric celestial bodies. New analytical results were 
determined to fill the gap between previous studies and astro- 
physical applications. In particular, we have derived generic for- 
mulas for the growth rate of the elliptical instability driven re- 
spectively by tides in non-synchronized bodies (TDEI) and li- 
bration in synchronized ones (LDEI), in the presence of imposed 
magnetic and thermal base fields. It was shown that an elliptical 
instability is strongly expected in the core of Europa, 55 CnC e, 
CoRoT-7b, and GJ 1214b, as well as in the subsurface ocean of 
Europa. Those results are valid for the present state of the con- 
sidered bodies and do not preclude any elliptical instabiUty in 
the past. For instance, the Early Earth core was clearly unstable, 
because of the larger gravitational distortions when the Moon 
was closer 

One can wonder about the signatures and consequences 
of such an instability on the planetary dynamics. A first con- 
sequence would be on the orbital evolution and synchroniza- 
tion process: indeed, the elliptical instability generates three- 
dimensional turbulent flows with cycles of growth, satura - 
tion, fluctuations, and relaminarization (e.g. lLe Bars et al.l2010h . 
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Fig. 8. Same as figure|6]but for CoroT-7b in figure (a), GJ 1214b 
in figure (b), and 55 CnC e in figure (c). 



Timescales involved range typically between the spin period and 
the growth time of the instability. Dissipation rates on the plane- 
tary scale, and consequently the orbital evolution, may then fol- 
low the same variations, with periods of rapid evolution when 
an elliptical instability is present, followed by more quiescent 
periods, for instance when the forbidden zone is reached. This 
increased dissipation s hould accelera t e the s ynchronization pro- 
cess, as described in iLe Bars et al.l (1201 Ol) . and this range of 
timescales should appear in the evolution of the spin rotation 
rate. 

The second consequence would be on heat flux variations a t 
the planetary surface. Inde ed, as shown in lCebron et al.l(l2010ct) : 
iLavorel & Le Bars! (I201 0'). flows driven by elliptical instability 
are very efficient in transporting heat by advection. As a result, 
subadiabatic cores should not be regarded as thermal blankets 
in the presence of the elliptical instability. The usual Nusselt 
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number Nu - QioilQdiff, the ratio between the total outward 
heat flux Q,oi and the purely diffusive outward heat flux Qdiff, 
which is associated to this heat advection, is given by the follow- 
ing sc ahng law, verified both experimentally (iLavorel & Le Bars! 
1201 Ol) and numerically JCebron et al.ll2010cl) : 



Nu 



0.01 



(37) 



This leads to a total outward heat flux advected by the ellipti- 
cal instability about Nu ^ 3-10"* times greater than a purely 
diffusive outward heat flux. Besides, in the presence of natural 
thermal convection, the superimposition of chaotic elliptically 
driven flows would induce large-scale variations in the same am- 
plitude. 

Finally, internal flows driven by elliptical instability are di- 
rectly responsible for magnetic field generation. The question 
of whether LDEI and TDEI are dynamos-capable is still open 
and remains out of the reach of the currently available numeri- 
cal capacity, but elliptically driven flows induce a magnetic field 
from an existing background one. To estimate a typical ampli- 
tude of such an induced field, we can use the results of our WKB 
approach. This shows that the dimensionless induced magnetic 
field inside the core and the instabiUty velocity uo are systemat- 
ically related by 



B = i 



Rm kr 



k^ 



Uo, 



(38) 



where k and k^^ are the norm and the axial component, respec- 
tively, of the wave vector of the excited mode of the ellipti- 
cal instability. This generic expression shows that the induced 
magnetic field due to the elliptical instability is systematically 
proportional to and in quadrature with the velocity field due 
to the instability. For the TDEI, fe^, - k/2, and for the LDEI, 
kx, = cJo k/4 (see respectively appendices lAl and iBli. Then, as- 
suming that at saturation, the typical flow induced by instability 
is commensurate with the differential rotation between the fluid 
and the elliptical distortion (see TablelTJ, we estimate the surface 
field by 



^ surf 



Rmkx^ aQ.(R2 



fe2 



Aii IR2V 

1i\r) 



(39) 



where R is the planet radius. Starting from the Jovian magnetic 
field component along the rotation axis, the LDEI in Europa sub- 
surface ocean is capable of inducing surface variations of up 
to ~ 0.1 % of the ambient field (reached for k - In /{I - 77) 
and optical librations, i.e. e - 0.0188), and LDEI in its core 
up to ~ 100 % of the ambient field at the surface (reached 
for k = 27r/(l - rf) and physical librations, i.e. e - 2 ■ 10'^). 
Cons i dering Galileo's E4 flyby of Europa (see IZimmer et al.l 
l2000t iKabin et aP 11999 ). the background z-component of the 
magnetic field is modified from 410 nT to 380nT at a distance of 
Europa about 1.5 Europa radius. This modification of 30 nT is 
on the same order of magnitude as the possibly LDEI induced 
magnetic field from the core. Internal s ources should thus be 
considered in addition to plasma currents (iKabin et al.lll999h for 
interpretating Europa's magnetic signal. 

To finish with, one should notice that the TDEI and LDEI 
studied here for telluric plan ets, can also affect the g iant gaseous 
planets of our solar system (IWicht & TilgnerluOlOl) . as well as 
extrasolar gaseous planets such as the hot-Jupiters, whose dra- 
matic tidal deformations should excite vigorous elliptical insta- 
bilities in the planetary atmospheres, but als o in th eir host stars 
(.Rieutord.2003: .Ou et al...2007; .Cebron et ani20Ill) . 
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Appendix A: Stability analysis for 
non-synchronized systems 

Solving ( fT6b . we find the trajectories 



Vl +yS COS(f yi -/3^ ) 

yfT^ sinit ^/^^) 




(A.l) 



where without loss of generality, the origin of time has been de- 
fined such that the initial position along the trajectory is X(,=o) = 
[xi,X2,X2] - [r sjl +/3 ,0,0] (the results of the WKB analysis 
do not depend on this chosen initial position). 

The solution to equation (l22T i at order 1 in yS along a stream- 
line (lA.lb writes as 



k = 



^.v. 



kx,,, cosf-fci-,„ smt+/3kx^^^ sinf 

kx,^ sint + kx^^ cos t+ 13 kx^,, sinf 

ko cos a 



(A.2) 



where k^^^, k^^g, and ko are constant, and a is the angle be- 
tween the wave vector and the rotation axis. We define the 
phase (p of the wave vector by writing kx^„ - kosma cos cp, 
kxio ~ ^osina sin (p. 

We are now in a position to solve the system of linearized 
equations given in section l33] To do so, we use as unknowns the 
vertical velocity of the perturbed field ut, and its vertical vorticity 
W3 = dxiU2-dx2U\ = i(kx^U2-kx2Ui), as well as the vertical com- 
ponent bi of the perturbed magnetic field and the corresponding 
magnetic vertical vorticity C3 = i(kx,b2 - kx^bi). The resolution 
is then straightforward (see 'Herreman et al.l l2009l) . 

At order in /3, the system reduces to an harmonic equa- 
tion for M3, giving a dispersion relation with a pulsation / = 
2 (1 -H Q'-') co s Q, with the qua n tity Q'-' = Q o,f,/(0„„n - Qgrh) al- 
ready usedbv lKerswelll(l2002l) : lLe Bars et a l. (2010). Solvability 
conditions imply non-trivial solutions only if / = 1, which gives 
the resonance condition cos a = 1/(2 (1 + QP)) € [-1,1]. This 
means that the instability cannot grow when Q./Q.orb G [- 1 ; 1 /3], 
which is the so-called forbidden zone. Outside this band, the 
growth rate is determined by the nullity of the determinant of 
the solvability condition system. It is then maximized over all 
values of wave vector phase 0. The maximum is obtained for 
(p = n/4 and the inviscid growth rate writes as 



(2QC + 3)2 



^Ii' 



'"" 16|l+fiGp 

k^ Ra r d^O 



4^2 



A 



8(1 +k* E^IPr^)\\+QP 



4|l+QC|3(i +Rm^k-^) 
(A.3) 



with 



f = 



(/fc4 + Rm^) Ra rdrO-lk^ RmA(\+k'^ E^/Pr^) 
(2QC + 3)2 (1 + ;^ E^/Pr^) (k^ + Rm^) ' 



(A.4) 



where r and dr@ are the radius and the dimensionless radial 
gradient of temperature base field on the considered streamline, 
respectively. In the absence of viscous boundary damping (dis- 
cussed in section 13.6b . the inviscid growth rate is a correct ap- 
proximation of the viscous growth rate when the viscous dif- 
fusion term -k^E in equation ( fT9l ) is negligible, i.e. for pertur- 
bations of wavelength greater than the Ekman thickness Ve. 
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Equation ( IA.3b takes the thermal and magnetic diffusions into 
account. In a typical liquid core, the thermal Prandtl number is 
about (9(0. 1-1). Neglecting the viscous diffusion term also leads 
to neglecting the thermal diffusion term -k'E/Pr of equation 
( fT9] l. But since the magnetic Prand tl number is about Pm x 10"^ 
(e.g. IChristensen & AubertI 120061) . the magnetic diffusion term 
-k^/Rm - -k^E/Pm of equation (ISTT l appears as the dominant 
diffusion mechanism. It can nevertheless be neglected within the 
limit of large wavenumbers, where equation (IA.3) gives equation 
(El. 



Appendix B: Stability analysis for synchronized 
systems 

We consider that the differential rotation between the fluid and 
the elliptical distortion writes as 1 - y^t) = ecos(wof + q). 
The phase q is introduced here because we define the origin of 
time such that the initial position of the considered trajectory is 
X2 = 0. Equation (fTSI l can be solved analytically to find the tra- 
jectories. In the particular case ^ = 0, the solution can be written 
in the following compact form: 



Vl +P COS — sin(w„f) Jl -B^ 
Vl -P sin — sin(wof) Jl -0- 



(B.l) 



with Sg = sgn(cos((yof))- In the general case, the wave vector 
associated to this flow writes as 



k ■ 
kx2o + -j^isin(aJot + q)-smq)(l+/3)e 

ko cos a 



(B.2) 



where k^^^, k^^^, and ko are constant and a is the angle between 
the wave vector and the rotation axis. We define again the phase 
(phy kx,g = kg sin a cos (p,kx-,„ - K sin a sin (p. 

At leading order in e/?, the dispersion relation gives / - 
2 cos fl, and the solvability conditions system admits non-trivial 
solutions for / - a»o/2. Consequently, the authorized band is 
given by cos a = (Ool^ G [-1, 1] i-e. |Wo| < 4. The growth rate 
is determined at order 1 in eyS and must be maximized above the 
phases q and (p. The maximum is reached for q -Q,(p - jt/4 and 
gives in the absence of thermal and magnetic fields 



Ojm. 16 + (ajI 



64 



(B.3) 



iKerswell & Malkua (Il998h have performed a global ap- 
proach to the same instability, explicitly considering inertial 
waves coupling in a spheroidal geometry. In the absence of 
magnetic and thermal fields, they found a maximum inviscid 
growth rate (T,„,,/(e P) - 25/128, very close to our value 
o'invli^ P) - 17/64. The small difference between the two val- 
ues is due to the influence of the spheroidal geometry consid- 
ered in IKerswell & Malkijsl (Il998h . leading to more restrictive 
conditions for de stabilization than our local analysis. Similarly, 
for purely hydrodynamic flow with a stationary deformation as 
studied in section 13741 our analysis gives the inviscid growth rate 
Cinv/P = 9/16, whereas a global analysis with inertial waves of 
a spheroid leads to the slightly lower value cr,„,,//3 =1/2 (see 



iLacaze et al.ll2004l) . Taking the uncertainties on the different pa- 
rameters for planetary application into account, these small dif- 
ferences can be disregarded and the local approach can be used 
confidently, which presents the strong advantage of providing an 
explicit formula for the growth rate. 

Taking a buoyancy of order e/3 into account, as well as the 
induction equation and a Lorentz force on the order of e/? in the 
presence of an imposed vertical magnetic field Bq, we obtain the 
growth rate 



16 H- w^ 



^inv ~ 



with 



6 = 



64 



^ ^(eP)^-4ajU2' 



coiA 



(16 + ajl)k^Rardr0 
16 (1 + k-'^colRm^lA) ~ 16 (w^ -h 4 fe^ £2/p^2) 



(B.4) 



(16 -I- w2) (4 ;^ -H w2 Rm^) (ojI -h 4 /fc^ E'^jPi^) 



(B.5) 



where F = (64 H- uli^k^ + Rm^(l6 + ul)))Ra r dfi - 
4-k^colRmA((juil + 4-k'^E^/Pp'), and where r and drO are respec- 
tively the radius and the dimensionless temperature radial gra- 
dient of the considered streamline. As discussed in appendix lAl 
in astrophysical applications, the thermal diffusion can be ne- 
glected (E/Pr - 0). In this case, equation ( IB.4I ) gives equation 
(|25] | in the limit of large wavenumbers. 



Appendix C: Is the diurnal tide stabilizing or 
destabilizing for the elliptical instability? 

In this work, the periodic forcing due to diurnal tides has been 
neglected for synchronized bodies. It is thus legitimate to calcu- 
late its influence on the LDEI growth. To answer this question, 
we consider the simplest but severe case of a body with a non ro- 
tating (i.e. y - 0) diurnal tide of amplitude /?(,) = /?i cos(Mf + q), 
where /5i - 3e/5 and there is no global rotation. Then, the base 
flow (fTTT l reduces to 

U = [-(1 +ySi cos(Mf + q)) X2 Cxi 

H-(l -ySl COs(Mf -H^)) JCl CxJ (C.l) 

Once again, the phase q is introduced here because we fix the 
phase of streamlines. The streamlines are not known analyti- 
cally. The dispersion relation gives the pulsation f - 2 cos a, 
and the solvability conditions give resonances for f - (2 + M)/2 
and / = (2 - M)/2. The authorized band is thus \M\ < 6. In 
the limit of low M, the maximal growth rate is obtained with 
/ = (2 - M)/2 for 4> = -7r/4 and ^ = 0: 



^inv ~ 



51 



— + — M+ -- 
16 64 



1 I5n^ 



8 32 



M' 



Pi 



(C.2) 



at the order O(A^) + 0{J3\ M^) + 0{fi\ ). This expression agrees 
with equati on (l24l l for A = and Ra - 0, in the Umit M -0. The 
expression (IC.2I I shows that slow oscillations of the amplitude of 
the tidal bulge are not inhibiting for the elliptical instability. On 
the contrary, the growth rate is enhanced compared to the case 
of constant amplitude for low M, which means that the diurnal 
tide would be destabilizing in this case. In the case of planetary 
interest with M - \, the maximum growth rate writes as 



25 
^'■""-l28^'' 



(C.3) 
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which again shows that the diurnal tide can drive an elhptical 
instabihty. This effect will thus be superimposed on the TDEI 
and LDEI mechanisms already studied, but with a slower growth 
rate, since ySi <K /?. 



Appendix D: Resonances of gravito-inertial waves 

To clarify the influence of a thermal field with a buoyancy force 
that is on the order of zero in /3, we consider the generic case with 
the elliptical gravitational iso-potentials of ellipticity n p and el- 
liptical base-field isotherms of ellipticity m /3, where n and m are 
arbitrary constants. This generic no tation is needed to deal wit h 
all cases, as s tudied for instance in 'Le Bars & Le D izisI (l2006h . 
iLavorel & Le Bars (201^) and [Cebron et al. (2010c'rin focus- 
ing on dynamic tides in a non-synchronized system, one would 
expect the isotherms to follow the streamlines (because of the 
small thermal diffusion coefficient), as well as the iso-potentials. 
Hence, n - m - 1 . In contrast, looking at the elliptical instability 
in a subsurface ocean underlain by a rigid mantle, one would ex- 
pect the iso-potentials to remain quasi-circular (n = 0). Besides 
this, in the presence of a static bulge, one would expect the sys- 
tem to naturally return to a configuration with m = by the gen- 
eration of baroclinic motions. All situations with < m, n < \ 
are possible, and one can even imagine other azimuthal period- 
icities, for instance those due to local variations of temperature. 
In the case of the TDEI, the WKB approach including 
zeroth-order buoyancy forces in /? is tractable. Unlike the cases 
studied in appendices |A] and |Bj the forbidden band where the 
instability does not exist is now modified by the thermal field 

and is given by /o < 1, where /o - V^jTTTy^pT'^flTs^. 
Outside this band, the inviscid growth rate is given, in the limit 
of large wavenumbers, by: 



{2QP -H 3)2 -h [1 + 2(1 + QP){m - n) - n] Ra rdr& 



yS 



16 |1 +QG|3+4/fQ|i +QG| ^g^Q 
A 1 



<D.l) 



4|1 -HQGp /3' 



We can compare the role of the temperature field in (1241 and 
( ID. 11 1. In expression (l24l i. the temperature field acts as a simple 
supplementary stabilizing term that corrects the inviscid purely 
hydrodynamic growth rate. But in the derivation of equation 
dD.lb . waves and resonances (as well as the forbidden band) 
are modified by buoyancy forces, leading to a modification of 
the prefactor of /3. Actually, the elliptical instabihty now results 
from resonances of gravito-inertial waves and should be called 
the g ravi to-elliptical instability (see Le Bars & Le DizesI 120061: 
iGuimbard et al. 2010) . As shown below, the supplementary res- 
onances associated to gravito-inertial waves allow the tempera- 
ture to be destabilizing in certain cases. The same conclusions 
can be obtained for the magnetic field when Lorentz forces are 
taken into account at zeroth order in /? in the limit of ideal magne- 
tohydrodynamic: the elliptical instabilit y then results from reso - 
nances between magneto-inertial waves (lKerswellll993al[l994l) : 
the forbidden band is modified by the magnetic field; and as 
shown by Lebovitz & Zweibel (2004); and the magnetic field 
can be either stabilizing or destabilizin g dep e nding on the case 
being considered (see also .HerremanI l2009t JMizerski & Baierl 
l2009l 1201 lb . Naturally, these conclusions are also valid for the 

LDEI.^ 

In Le Bars & Le DizesI (l2006h and in the experiments of 
ILavorel & Le Bars! (1201 Ol) . the TDEI is studied for a station- 
ary bulge (Q'^ = 0) with circular iso-potentials and elliptical 
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Fig.D.l. Growth rate of the TDEI for a cylindrical shell of 
aspect ratio H/R = 2 with an elliptical cross section (tj = 
0.2, E = 0.0036, yS = 0.47, Pr = 1) and an ellipsoidal shell 
(77 = 0.3, E = 0.0029, /3 = 0.317, Pr = 1) with a rotation 
axis of length c - (a + b)/2. The figure compares the numeri- 
cal growth rate in the autogravitating case where the gravity is 
given by the Poisson eq uation for the gravita tional potential of a 
homogeneous fluid (see lCebron et al.l2010d, for details), and the 
case where the gravity is played by a centrifugal force (« = 0), 
as in the experiments of Lavorel & Le Bars (2010). The numer- 
ical growth rate is translated vertically in the figure, to match 
the inviscid growth rate 9/16 at Ra - 0, which corresponds to 
a surfacic damping term coefficient a - 3.24 for the cylindrical 
shell (squares and triangles) and a = 3.1 for the ellipsoidal shell 
(circles and diamonds). Theoretical growth rates are shown by a 
continuous red line for n = 0, m = 1 and a dashed blue fine for 
n — tn — 1 . 



isotherms (n - 0, m = 1), and the considered temperature pro- 
file gives r d,® = - 1 . In this particular case, in the absence of 
magnetic field, equation (ID. lb recovers their result: 



9-3Ra 

^/3. 

16 - 4 Ra 



(D.2) 



As already noticed by Le Ba rs & Le D izes (2006), a thermal sta- 
ble stratification (Ra < 0) is then destabilizing for the elliptical 
instability, but, in constrast, the temperature field stabilizes the 
instability for n = m = 1 . This high sensitivity of the growth 
rate of the elliptical instability to the specific gravitational and 
thermal fields is con firmed by numerical simulations. Using the 
method described in lCebron et al.l (1201 Od) . we consider the sim- 
ple case Q*^ - and K - 0, the temperature field being es- 
tablished by a temperature contrast between the two boundaries. 
When n - m - \, the growth rate (ID. lb is enhanced when Ra 
is increased. As shown in figure ID. II this is in perfect agree- 
ment with the numeric al simulations in a cylindri cal shell. In the 
experimental setup of ILavorel & Le BarsI (l2010l). the gravity is 
replaced by the centrifugal acceleration, as in'Carrig an & Bussa 
(1983), and the associated equipotentials are circular, i.e. n = 
and m = 0, as shown in figure IdTI In this case an increasing Ra 
indeed leads to a lower growth rate and the numerical simula- 
tions agree with the predicted growth rate. 

The conclusions for an autogravitating ellipsoidal shell are 
more complex. Cebron et al. (2010c) find that an increasing Ra 
leads to a lower growth rate, which contradicts the prediction 
of the theoretical growth rate. This difference comes from the 
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fact that the thermal stratification propagates the influence of the 
boundary inside the bulk: the WKB analysis, based on local sta- 
bility, cannot handle this feature. In the spherical geometry, we 
can, however, notice that changes induced by Ra are small for 
-I < Ra < 1 and remain close to the estimates for an autogravi- 
tating cylinder. 

Those results clearly illustrate the high sensitivity of the 
growth rate of the elliptical instability to the specific gravita- 
tional and thermal fields, as well as to the considered geometry. 
In planetary applications, stratification (i.e. Ra < 0) generally 
leads to stabilization, as in the limit of small Ra presented in the 
main text; however, in this case, stratification can only stabilize 
elliptical instability when Ra - (9(1 - 10), which is never the 
case. 
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